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Abstract 

We present a new program package for calculating one-loop Feynman 
integrals, based on a new method avoiding Feynman parametrization and 
the contraction due to Passarino and Veltman. The package is calculating 
one-, two- and three-point functions both algebraically and numerically 
to all tensor cases. This program is written as a package for Maple. An 
additional Mathematica version is planned later. 

1 Introduction 

In particle physics the calculation of one- and higher loop corrections to parti- 
cle processes is necessary to keep track with the increasing accuracy of particle 
colliders. In this context it is very useful if a fast and automatized calculation 
of the contributing Feynman diagramms is provided by computer. On the most 
basic loop correction order, the one-loop level, several attempts for such pro- 
grams are already available, like FF [Jl3| - calculating scalar one-loop integrals 
in Fortran - or FeynCalc Jl2| - a Mathematica jll| package - which exhibits 
the usual Passarino- Veltman coefficients of tensor integrals. Both programs are 
using the standard procedure for calculating one-loop diagrams, namely Feyn- 
man parametrization, Wick rotation and the Passarino- Veltman contraction for 
tensor structure (cf. || and ||). 

Although this standard procedure became very popular during the last 
decade, it turned out to be very inconvenient for being implemented on computer 
if the integral involves higher tensor structure. For that reason we investigated 
in a new method for calculating arbitrary one-loop integrals which is already 
published for the two- and three-point case (cf. jl| and ||). 

Our approach is completely different from the standard procedure and the 
structure of other packages. We are not forced to invert big matrices which is 
the most difficult step in computing the tensor structure of the integrals. In 
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our method the loop momentum integrations are performed directly, using the 
residue theorem and introducing 1Z functions (see ||), a class of special functions 
with some nice and useful features. Instead of big matrices we get some recursion 
relations for the involved 1Z functions (cf. Q or [Q), which allow a faster and 
more memory efficient evaluation of the integrals under consideration. 

The program is designed to reduce the TZ functions to a set of fundamental 
functions which are rewritten in terms of logarithms and dilogarithms. Up 
to now the different one-, two- and three-point functions are programmed as 
procedures in Maple fll(i|| . Therefore our results are available numerically as 
well as algebraically. We tested our procedures using Maple V on a VAX 4000- 
90 Workstation. 



2 Description of the method 

A detailed description of the method may be found in (!]] or Q . Here we would 
like to give a brief introductory overview only. The calculation is based on 
the well-known dimensional regularization method, as described for example 
in j(| . Instead of introducing Feynman parameters the integral is split up into 
orthogonal and parallel space (cf . || ) , so that for example the two-point integral 
takes on the form 

B iPoPl) (q,m u m 2 ) = 



r(^i) J "7 x V1P2 







with 



Vi = (Z|, +qf - ll-mf + ig 
~ L \\ 



V 2 = lf\-l±- ml + ig 



Here we present the most general case, depending on one external momentum 
q and two arbitrary masses m\ and m,2- The splitting of space requires a dis- 
tinction between In and Zj_ which we define as 

I 



'11 = -?= ; l ^ = \ l l- 12 ( 2 ) 



In this notation 2 11 describes the component of I which is parallel to the external 
momentum q, whereas l± represents the orthogonal complement. 



The integrals (|lj) can be reduced to a sum of 1Z functions by using ( |A.l[ ) 
with the help of partial fraction and residue theorem. The complete procedure 
is described in detail in [0 and 0]. The result in the two-point-case may be 
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written according to (^]: 



EI 



B ( - PoPl \q 1 m u m 2 ) = in 



3 e (-1) 2 ( e — )- 



K 



q{-2e+ Pl + 1) 

r(|- e) 



Pi+boi {-H 1 + e > i; -™? + iQ, -(§ + A) 2 ) 



Po 



+ £( , ?)a 9 fl( e -*l±M,M + i)x 



2 2 

^_ e+ £i±w (-^ + e, 1; -m| + i e , -(f - A) 2 ) 



(3) 



where e = 4 2 " D represents the usual dimensional regularization parameter. We 
used the abbreviations 

ai — (_l^Pa-[pa] _|_ ^1+po-bo] 



.4 



9 9 
777-2 _ m l 



2 9 



i, z even 

i + 1, z odd 



The result of the three-point integral is obtained by a quite similar procedure 
and may be found in M or M . 



3 Evaluation and expansion of the 1Z functions 

After having expressed all possible integrals in terms of 1Z functions we get easy 
and fast computable results by applying some recurrence relations for the 1Z 
functions under consideration. 

We achieve these relations by first reducing the number of different 1Z func- 
tions to a set of fundamental 1Z functions. This is possible due to the already 



mentioned features of the 1Z functions, especially the relation ( |A.2| ) and (A. 3). 
The reduction takes place in two steps. First of all we increase the parameter 
of the 1Z function to get it close to or 1. This can easily be achieved by us- 



ing formula (JA. 2|) . In the second step we reduce the index of the 1Z function 



by virtue of a convenient formula derived from (A. 3). In the case of the two 
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point-function it has the form: 



7^(61,62; Zl,Z 2 ) 



-1 



t + /3-l 
[ {{t - l)zxz 2 ) U t - 2 (bi,b 2 ; z 1 ,z 2 ) 

+ ((l-t-b 1 )z 1 + (l-t-b 2 )z 2 ) 
x Kt-i(bi,b 2 ;zi,Z2)} 



(4) 



Now the power of this method is obvious. Both steps mentioned above may be 
performed easily by computer, due to the recursive definition of the formulae 
( |Al|) and flOl ). 

Up to now we have acquired a result which is exact to all powers of e, 
though we are interested only in an expansion in powers of e. Therefore we 
now substitute the expansion for the fundamental TZ functions in terms of e, 
resulting with simply computable functions like logarithms and dilogarithms. 
This is a rather cumbersome task but has to be performed only once for any 
N-point function. The procedure, especially its subtleties concerning analytical 
continuation, is described in detail in another paper [Q. For completeness we 
here just cite the expansion for the scalar two-point case 



Tl- e (-h + £, 1, zi, z 2 ) 



iri (z 2 -z 1 )i 



B{e,\) 



1 -s / Z l 
+ 4;Z- 



"2 "1 



z 2 



1 + 4/1- 



-l+2e 



1-W1-^ 

z 2 



-l+2e' 



+ 0{e 2 ) 



The three-point integrals 

C (MPlp2) (gi,«2,m 1 ,m 2 ,m 3 ) = x 

1 \ 2 ) 

OO OO OO 



(5) 



(6) 







Vi V 2 Vz 



V x = {k\\+<IxT -l\-l\-m\ + iQ\ 
V 2 = (l Q \\ + q 20 f - (h\\ + q 2 if m\ + ig 

^3 = - -l±-ml + ig 

are calculated in a quite similar way (cf. or 0]). Again we treat the 
general case with two independent external momenta q\ , q 2 and three masses 
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mi, TO2, TO3. For the momentum components we write in a quite natural gener- 
alization: 



foil - 



'9 2 



, 91 ■ 92 

9 2 = 92 — 9l 

9i 



920 



2 - / 2 

L|| ' 

91 ■ 92 



921 = \/gfo -9l 



The evaluation is resulting in the fundamental 7?. functions 

ft_ 2e (e, £,!;&, 2/, z) = 



{ 



1 



2e z 



Li 2 (l - ~) + A(x,z) 
+Li a (l-|) +A(y,z)]} + 0(e 3 ) 
7e_i_2 £ (e, e, 1; 1, 2/, z) = (1 - 2e)z 1 - 2e 
+ e(x + y) + 2e 2 — a; In a; — ylny + zLi 2 ^1 J 

+z A(x, z) + z Li 2 - -) + z A(y, z) + 0(e 3 ) 



here we introduced the A-Function. 
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Table 1: Input and Output for the one- 


-loop procedures 




input 


output 




OneLooplPt(p, to) 

OneLoop2Pt(po,Pi, 9, mi, TO2) 

OneLoop3Pt (p , Pi , P2 , qi , 920, 921 , mi , ms , ms) 


[e~ 1 -term, e° 
{e~ 1 -term, e° 
[s~ 2 -term, e~ 


-term] 
-term] 

1 -term, e°-term, {abb.}] 



A(x, z) = In ( 1 ) 77 ( x, 



- In ; 



X] I x — z, 



-ri\x-z, — 

z 



(9) 



Note that our relations are valid in all kinematical regions, so that the results 
of apply to all the cases of our program. For a detailed discussion sec j|J . 
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4 The One-loop functions 

Our aim is to calculate one-loop integrals both algebraically and numerically. 
Therefore the algorithm is implemented in Maple a 'language' for symbolic 
computations. To make the implementation for any user most comfortable, the 
procedures are bound to a package, which may be read from the Maple system 
during run-time. The user-interface consists of three elementary functions 

• OneLooplPt(p, m) 

• OneLoop2Pt(poiPii 9> m i, m 2) 

• OneLoop3Pt(p ,Pi,P2, 9i, 920, <?2i, m-i, m 2 , m 3 ) 

The functions directly correspond to the integrals For completness the 

one-point function 



is implemented. 

We followed again our notation which distinguishes between parallel and 
orthogonal space. The masses are denoted by nij and the definitions of (||) and 
(0) for the momenta are of course still valid. We would like to emphasize the 
fact that in our notation the indices Pi represent the powers of the different 
components of the loop momentum I and should not be mixed with Lorentz 
indices of the corresponding tensor. This notation is perhaps not very often used, 
but it corresponds directly to our method of integration. The usual notation 
due to ||] may be recovered by an additional algorithm also cointained in our 
package. 

The results of the one-loop functions are Laurent expansions in terms of the 
ultraviolet regulator e. The output of the described procedures consists of a list 
where the significant coefficients of this expansion are denoted (see tabic . 

In the output "e _1 -term" is meant to represent the divergent part (the co- 
efficient of i) whereas "e°-term" is describing the finite part. Of course the 
£ -2 -term has to vanish in any case, but is left in order to check the accuracy of 
the program in numerical calculations. In the three-point case "abb." is written 
for some abbreviations which are introduced to get a more compact result. 

As long as the arguments of the OneLoop functions are symbols the output 
is given algebraically, whereas of course numbers inserted for the arguments 
imply a numerical result. Moreover, in the numerical case, the program sets the 
value of g, the imaginary part of the propagators, five digits higher than the 
numerical accuracy of the whole calculation. In the algebraic case just 'rho' is 
returned. 

For different calculations of the same OneLoop function it is of course rather 
inconvenient to start the whole program again. Therefore our package includes 
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the functions OneLoopLiblPt(n), OneLoopLib2Pt(n) and OneLoopLib3Pt(n) 
which generate a library of all OneLooplPt, OncLoop2Pt and OneLoop3Pt 
functions respectively up to the tensor rank n. Each function is written to one 
file. The filenames correspond to the integrals in the following way: 



TWO00.RES OneLoop2Pt(0, 0, q, mi, ra 2 ) 
TWO10.RES OncLoop2Pt(l, 0, q, m u ra 2 ) 
THREE000.RES OneLoop3Pt(0, 0, 0, qi, q 20 , 

q 21 ,m 1 ,m 2 ,m 3 ) 
THREE100.RES OncLoop3Pt(l, 0, 0, q u q 20 , 

q 2 i, mi, m 2 ,m 3 ) 
THREE010.RES OneLoop3Pt(0, 1, 0, q u q 20 , 

q 2 i, mi, m 2 ,m 3 ) 



The path of the library may be specified with the variable "tensorpath" . 



5 An example 



To demonstrate how our package works we illustrate the calculation of the 
electron self-energy in QED. Starting with QED Fcynman rules we get: 



ONE0.RES 



OneLooplPt(0,m) 





2 



/ 



d D l 7M (^ + ; + m ) 7 M 

(2ir) D [{l + q) 2 -m 2 }l 2 




J (2n) D [(l + q) 2 -m 2 ]P 
. o f d D l 1 
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x \l vl v {2 -D) + q„j u (2 -D) + mD] 

Now we split into orthogonal and parallel space noting that the 7 matrix vector 
is split into 7|| and 7^: 

= -j^yj{(2~D)[ lv B^(q,m,0) 

- 1± -B^\q,m,0)} 

+ (2-D)q- lll B^°\q,m,0) 

+DmB^(q,m,0)} 

At this stage of calculation we now take the opportunity to incorporate our 
procedures which solve the ocurring integrals B^ Po ' Pl \ To be specific we have to 
compute the functions OneLoop2Pt(l, 0, q, m, 0), OneLoop2Pt(0, 1, q, m, 0) and 
OncLoop2Pt(0, 0, q, m, 0). If we restrict ourselves in this small exercise to the 
divergent part C(e _1 ), we end with the following results: 

OneLoop2Pt(l, 0, q, m, 0) = £ (10) (g, m, 0) = 

V {"'} 

OncLoop2Pt(0, 1, q, m, 0) = S (01) (q, m, 0) = 
OneLoop2Pt(0, 0, q, m, 0) = B m (q, m, 0) = 

■ 9 2 

So our final result is 

e 2 

= (-9 ■ 7|| + 4m) + finite 

e 2 

= T7T^(-^ + 4m) + finite (11) 
which is the well-known result for the divergent part of the electron self-energy. 



2e Q ' 4q 3 



6 Conclusion and further projects 

Our aim was to introduce a program package that gives any computer user 
the possibilty to calculate one-loop particle processes without much effort in 
manpower and hardware. 

As a further project we plan to extend this package to the four- and five-point 
function. Moreover we want to export our program package to Mathcmatica 
[ p"T| , because from our point of view Mathematica is the more common language 
than Maple. Planned is also a menu-driven user interface. But this might be 
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rather difficult, because we haven't yet found a convenient C-language package 
for algebraic calculations. (For hints contact L. Briichcr) 



Appendix 

Conventions 

D = 4 — 2e Dimension in the dimensinal 

-regularisation-mcthod 
TZt(bi,b2, zi,z 2 ) TZ function, t is called index, 

61,62 parameters, 

z\ , z 2 arguments 
(3 = J27=o hi sum 01 parameters of an 

71 function 

Some fundamental properties of 7Z functions 

• formula for integrals 

k 



[ x"- 1 f[(zi + Wi x)- bi dx 

1=1 

k 

= B((3-a,a)7l a . (b,-)l[w^ (A.l) 



1=1 



increasing of parameters 

- tz t K t -i(b + e t ,z) (A.2) 



associated functions (k is the number of parameters) 

k 

£ 4 (fc) (^b,z)7^ t _ 4 (b,z) = (A.3) 



i=0 

with 



t(/3-t-k) 



k 



:{t + i-Y,b jZj ^)Et\z) (A.4) 



REFERENCES 



10 



E\ (z) is the elementary symmetric function, gained from 

k 

{l + xz 1 )---{l + xz k ) = Y J X i E\ k \z) (A.5) 

i=0 
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